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The nature of source rocks of basaltic magmas plays a fundamental role in understanding the composition, 
structure and evolution of the solid earth. However, identification of source lithology of basalts remains 
uncertainty. Using a parameterization of multi-decadal melting experiments on a variety of peridotite and 
pyroxenite, we show here that a parameter called FC3MS value (FeO/CaO-3*MgO/Si0 2 , all in wt%) can 
identify most pyroxenite-derived basalts. The continental oceanic island basalt-like volcanic rocks 
(MgO>7.5%) (C-OIB) in eastern China and Mongolia are too high in the FC3MS value to be derived from 
peridotite source. The majority of the C-OIB in phase diagrams are equilibrium with garnet and 
clinopyroxene, indicating that garnet pyroxenite is the dominant source lithology. Our results demonstrate 
that many reputed evolved low magnesian C-OIBs in fact represent primary pyroxenite melts, suggesting 
that many previous geological and petrological interpretations of basalts based on the single peridotite 
model need to be reconsidered. 

For years it was commonly assumed that basaltic magmas are derived from mantle peridotite 1,2 . Although 
several lines of evidence from multi-disciplines indicate that pyroxenite (unless otherwise stated, pyroxenite, 
including amphibolite, refers to mafic lithology that have a similar mineralogy as peridotite, but contains less 
than 40% olivine. From compositional perspective, pyroxenite and peridotite generally have <30 wt% and >30 
wt% MgO, respectively, see discussions in Supplementary Information) may play a critical role in the generation 
of basaltic magmas 3 9 , there is uncertainty about identification of source lithology of basalts because the current 
criteria such as high Ni content, high Fe/Mn 10,11 and Zn/Fe ratios 12 , and low CaO content 13 in either melts or 
olivine, which might also be interpreted as a result of peridotite compositional diversity and variations of melting 
and crystallizing conditions of basaltic magmas 1418 (Supplementary Information). On the other hand, isotope 
and trace elements cannot determine whether basalts are pyroxenite-derived melts or the pyroxenite signature 
just imprinted onto peridotite, although they might tell us the important role of pyroxenite 1 ' -20 . Furthermore, the 
phase relations of basaltic melts, although can thermodynamically explain their overall compositional character- 
istics and melting pressures based on assumed source mineralogy, cannot constrain their source lithology and 
bulk composition 17,19,21 . 

Results 

FC3MS value can identify most pyroxenite-derived basalts. A large body of melting experiments on a variety of 
peridotite and pyroxenite at various melting conditions reported over the past three decades are parameterized in 
our study (Supplementary Figs. 1-5). In the commonly used simple plots, the pyroxenite melts, except for some 
low MgO (<5 wt%), high S0 2 (>60 wt%), high Ti0 2 (>6 wt%) and high FeO (> 15 wt%) melts, can overlap with 
peridotite melts if assuming that all the melts can be affected by olivine fractionation and accumulation. 
Importantly, the C-OIB can be grouped into either peridotite melts or pyroxenite melts in the simple plots 
and projections(Supplementary Figs. 1-4). This difficulty has also been realised by the recent review work of 
Lambart and co-workers 22 on the markers of the pyroxenite contribution in the major element compositions of 
oceanic basalts, and which would commonly make us believe that the pyroxenite and peridotite melts are mostly 
indistinguishable and identification of source lithology of the C-OIB seems impossible based on major elements 
(Supplementary Information). Surprisingly, however, most pyroxenite melts can be identified because of their 
higher FeO/CaO ratio (unless otherwise stated, FeO refers to total FeO) at a given MgO/Si0 2 ratio, and the two 
ratios in general showing positive and negative correlations for peridotite and pyroxenite melts, respectively 
(Fig. la). Almost all the peridotite and pyroxenite melts, at a given pressure, showing negative correlations 
between MgO/Si0 2 and FeO/CaO before clinopyroxene is melted out (Fig. 1). When volatiles such as H 2 0, 
C0 2 , CHO-fluid, and metasomatic minerals such as amphibole and phlogopite are present in peridotite source, 



SCIENTIFIC REPORTS | 3 : 1 856 | DOI: 1 0. 1 038/srep01 856 



1 




0.0 0.2 0.4 0.6 0.8 1.0 

MgO/SiO, 




o.o 1 ' — 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 
MgO/Si0 2 



Figure 1 | FeO/CaO versus MgO/Si0 2 of experimental peridotite- and 
pyroxenite-derived melts and the C-OIB. (a), The four red line is the 
upper limit of near-solidus melts compositions for iron-rich olivine 
websterite (HK-66), iron-rich lherzolite (PHN1611), and normal 
lherzolite (KLB-1 and KR4003), and of melts that equilibrium with 
orthopyroxene and olivine, respectively. The numbered blue circles 
represent near-solidus melts compositions of the olivine websterite and 
peridotites from 1 Gpa to 7 Gpa, and the black arrows refer to melting 
trends at constant pressure (Supplementary Figs. 8-10). The region 
depicted by red curve are volatile-free experimental peridotite melts 
with Mg # = 59.5—68.0. (b), The grey arrow shows the overall melting 
trend of pyroxenite (Supplementary Fig. 11). The arrows with white 
circles and yellow squares, olivine (Ol) and clinopyroxene (Cpx) 
fractination trends (10 wt% step) of two hypothetical parental melts 
representing primary alkaline and tholeiitic C-OIB (Supplementary 
Table 1). Olivne fractionation trends were obtained using PRIMELT2 
software 35 and assuming Fe 2 0 3 /Ti0 2 = 1 in that the C-OIB might have 
high oxygen fugacity. It is commonly difficult to evaluate clinopyroxene 
fractionation, here we assume equilibrium crystallization and the 
partition coefficients of all major elements are estimated from our 
collected experimental results (Supplementary Table 1). The region 
depicted by the black curve shows EPR-MORB compositions (7.5 wt% 
< MgO < 1 1 wt%) acquired from EarthChem database and only those 
total iron content reported data are selected. The grey areas show the 
global average compositions of MORB including original compositions 
(higher FC3MS value) and corrected primary compositions (lower 
FC3MS value) 23 . The C-OIB data acquired from GEOROCK database 
accessed December 2012 (Supplementary Table 2). Other data source 
given in Supplementary Information. 



the resultant melts FC3MS values are either equal to or lower than 
that of volatile-free peridotite melts (Fig. 1 and Supplementary Figs. 
12 and 13). The parameterized results of experimental peridotite 
melts shown in Fig. 1 also applied to theoretical calculations 
(Supplementary Fig. 13). In addition, the MgO/Si0 2 and FeO/CaO 
ratios of the peridotite melts in general are positively correlated with 
pressure (Supplementary Fig. 14), which is very similar as the 
pressure effect on FeO/MnO ratio (Supplementary Fig. 15), but the 
FC3MS value (FeO/CaO-3*MgO/Si0 2 , all in wt%) is pressure- 
independent (Supplementary Fig. 16). 

The near-solidus melts of iron-rich lherzolite PHN1611 is esti- 
mated with FC3MS —0.60 (Fig. 1), and which is found to be the 
highest FC3MS value for all the experimental Lherzolite melts 
(Figs. 1-2). The FC3MS value of normal mantle lherzolite melts is 
estimated less than 0.5 (Fig. 1). It should be noted that the FC3MS 
value of lherzolite melts is largely dependent on total alkali content 
and melting degree (Fig. 2a, c). The FC3MS values of most pyroxenite 
melts can significantly higher than peridotite melts at the same total 
alkali content, melting degree and Mg* (MgO/(MgO + FeO)*100, 
molar ratio), and can also higher than the near-solidus peridotite 
melts even if pyroxenite have very high degree partial melting 
(~20%-50%) (Fig. 2c). Noteworthy is the wide range of Mg* for both 
peridotite and pyroxenite melts, although it is commonly assumed 
that mantle-derived primary magmas have Mg* ~ 70-73. Never- 
theless, pyroxenite melts can have very broad range of Mg*, at least 
ranging from 40 to 70 (Fig. 2b). Therefore mantle-derived primary 
magmas would not necessarily have high Mg*, and which of course 
can show large variations in both Mg* and MgO content if we do not 
assume that basaltic magmas are equilibrium with typical mantle 
peridotite. Using the systematic variations of the FC3MS values ver- 
sus Na 2 0 + K 2 0 (wt%), Mg* (MgO/(MgO + FeO)* 100, molar ratio) 
and melting degree (F%) of a variety of peridotite melts (Supple- 
mentary Figs. 17-23), we determined the upper limit of the 
FC3MS value for olivine-controlled fractionated peridotite melts 
after —10-30 wt% olivine fractionation when MgO > 7.5 wt% or 
Mg* > 50 (Fig. 2b and Supplementary Fig. 20). Simply, for peridotite 
melts the upper limit of the FC3MS value is 0.65. Our statistical 
results show that the experimental peridotite- and pyroxenite- 
derived melts are -0.07 ± 0.51 (25, n = 656) and 0.46 ± 0.96 (25, 
n = 494) in the average FC3MS value, respectively (Supplementary 
Fig. 24a, b and Supplementary Table 3). 

Peridotite melts cannot explain the FC3MS value of most C-OIB. 

Interestingly, the C-OIB (unless otherwise stated, the C-OIB refers to 
samples that have MgO>7.5 wt% throughout the remainder of the 
paper) FC3MS = 0.68 ± 0.34 (25, n = 525), the alkaline C-OIB 
FC3MS = 0.65 ± 0.35(25, n = 356) and the tholeiitic C-OIB 
FC3MS = 0.70 ± 0.30 (25, n = 138) (Supplementary Figs. 24c, d 
and 25 and Supplementary Table 3), of which more than half samples 
have higher FC3MS value than the upper limit of peridotite melts, 
but all the C-OIB FC3MS values fall in the range of pyroxenite melts 
(Figs. 1-3). The FC3MS values of the C-OIB are not correlated with 
major element contents of Si0 2 , Ti0 2 , A1 2 0 3 and FeO, but showing 
weak negative correlations with MgO and CaO contents. In 
particular, the low Si0 2 (<43 wt%) and high FeO (>13.5 wt%) C- 
OIB show relative lower FC3MS values (generally less than 0.65), 
although having the most significant differences relative to peridotite 
melts in simple plots (Supplementary Figs. 1-3). Nevertheless, these 
samples still significantly differ from all the known peridotite melts 
when presented in Fig. 3a. Because there is an upper limit for the 
FC3MS value of the peridotite melts (Figs. 1-3), isobaric melting, 
polybaric melting, or melt mixing from various kinds of conditions, 
the resultant FC3MS values of peridotite melts still cannot higher 
than 0.65 (Supplementary Figs. 21-22). The high FC3MS values 
of these relative primitive C-OIBs might be explained by 
substantial amount fractionation of clinopyroxene or pyroxenite 
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Figure 2 | Diagrams showing the FC3MS value as functions of Na 2 0 + 
K 2 0, Mg* and melting degree (%). For clarity, only the C-OIB data and 
melting trends of thirteen representative peridotite and pyroxenite are 
shown. Detailed data point and data source given in Supplementary Figs. 
17-19. Grey areas in (a), (b) and (c) are peridotite melts (Mg*>60). The 
error bars shown in (a) and (c) are the standard deviation for average 
pyroxenite melts and normal and enriched peridotite melts. 

partial melting. Before presenting our preferred interpretation, we 
discuss why clinopyroxene fractionation of peridotite melts cannot 
reasonably explain the high FC3MS values of the C-OIB. 

For the peridotite model of basalts, it is widely accepted that alkal- 
ine basalts are lower degree partial melts than tholeiitic basalts. 
Here we assume that the —4% and —13% partial melts of iron-rich 



peridotite PHN1611 represent the primary alkaline and tholeiitic 
basalts of the C-OIB (Fig. 2c), and their FC3MS values are estimated 
~0.4 and —0.2, respectively (detailed hypothetical composition see 
Supplementary Table 1). For peridotite source the hypothetical 
FC3MS values can be regarded as the upper limit of the primary 
C-OIB, which higher than most primary MORB 23 (Fig. lb). In order 
to explain the FC3MS values and major element contents of the C- 
OIB, we need assume —0-15 wt% olivine fractionation combined 
with —30-40 wt% and —60-70 wt% clinopyroxene fractionation for 
the alkaline C-OIB and tholeiitic C-OIB, respectively (Figs. 1-3). 
Nevertheless, if the C-OIB are indeed evolved magmas the fraction- 
ating mineral should be mainly olivine rather than clinopyroxene. 
First, ascent and decompression of peridotite melts would cause 
olivine oversaturated and pyroxene undersaturated. It is difficult 
for peridotite melts with MgO>7.5 wt% to crystallize clinopyroxene. 
Many previous studies suggested that the C-OIB mostly contain only 
olivine phenocrysts or aphanitic 2427 , which are the characteristic 
features even for the C-OIB with MgO < 7.5 wt% 27 . In particular, 
the alkaline C-OIB generally carry a variety of mantle-derived xeno- 
liths, indicating that magma evolution is very limited during ascent 
in crustal conduits. Second, olivine is commonly the dominant frac- 
tionating mineral for basaltic magmas, including MORB and OIB, 
when MgO>7.5 wt% 13 - 28 . The C-OIB generally have higher alkali and 
lower CaO contents than MORB at a given MgO content (Supple- 
mentary Fig. 4), which would expand the stability field of olivine at 
the expense of clinopyroxene, giving rise to clinopyroxene crystal- 
liation in the C-OIB at lower MgO content than that for MORB. 

Although the high Fe/Mn ratios of the C-OIB might be explained 
by high pressure partial melting of peridotite or as characteristics of 
iron-rich or high Fe/Mn peridotite melts 1516 , the negative correlation 
between Fe/Mn ratios and Yb contents of the C-OIB in north China 
craton is inconsistent with significant fractionation of clinopyrox- 
ene 26 . Importantly, even for those C-OIB with MgO>10 wt% and 
those with Mg* > 60, their average FC3MS values are 0.54 ± 0.3 1 (28, 
n = 135) and 0.57 ± 0.32 (28, n = 238), respectively (Supplementary 
Table 3). Nearly half of these primitive C-OIB still have higher 
FC3MS values than the peridotite melts. Note that the wide range 
in the FC3MS value of the peidotite melts include Mg* = 60—68 
melts (Fig. 2b and Supplementary Fig. 20). Therefore, it is highly 
unrealistic to assume —30-70 wt% fractionation of clinopyroxene 
for explaining the high FC3MS values of the C-OIB. It should also be 
noted, however, that many previous studies using MgO > 6 wt% C- 
OIB to acquire primary melts composition that equilibrium with 
peridotite using olivine addition method 25,29 . Although the calcula- 
tion is not carried out in this study, the C-OIB with MgO > 6 wt% are 
significantly deviated from the olivine fractionation trends of peri- 
dotite melts (Fig. 2b), and even for the alkaline C-OIB the elevated 
FC3MS values require —50-70 wt% clinopyroxene fractionation, 
which not only rather unrealistic but also contradicts the initial 
assumption of olivine-controlled fractionation. 

The global MORB glass compositions (MgO > 7 wt%) show that 
the highest increase of the FC3MS values is 0.47 relative to corrected 
primary melts that equilibrium with typical mantle peridotite, and 
for MORB with MgO>7.5 wt%, the increase of FC3MS value is 
universally lower than 0.3 23 . Therefore if assuming that these 
MORB glass and the C-OIB are peridotite-derived and underwent 
a similar fractionation process, the tholeiitic (assuming primary melt 
with FC3MS = 0.2) and alkaline (assuming primary melt with 
FC3MS = 0.4) C-OIB then would have FC3MS values less than 0.5 
and 0.7, respectively. In addition, the EPR-MORB (MgO > 7. 5 wt%) 
wifhFC3MS = 0.33 ± 0.20 (28, n = 1510) (Supplementary Table 3), 
which is significantly lower than the C-OIB (0.68 ± 0.34 (28, n = 
525)). Noteworthy is the FC3MS values of the EPR-MORB increase 
from —0.35 to —0.65 when MgO content decrease from 8.5 wt% to 
7.5 wt%, which accompanied by decrease of A1 2 0 3 and CaO contents 
and increase of FeO and Ti0 2 contents. These variations are very 
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Figure 3 | FC3MS values versus major elements contents and La/Yb of the C-OIB and experimental peridotite- and pyroxenite-derived melts, (a), The 
blue and yellow arrows showing the fractionation (increasing the FC3MS value) and accumulation (decreasing the FC3MS value) trends of olivine. The 
blue and yellow filled squares showing the fractionation trends of clinopyroxene, of which the larger squares represent the hypothetical primary magmas, 
blue, alkaline basalt, yellow, tholeiitic basalt. The brown filled circles, EPR-MORB compositions (7.5 wt% < MgO < 1 1 wt%). The red filled diamond, 
C0 2 and CHO fluid present peridotite melts, (b), Group 1, Fujian tholeiitic basalt™. Group 2, Tariat alkaline basalt 27 . Group 3 and group 4, Shuangliao 
alkaline basalt 31 . All the four groups data see supplementary Table 4. Group 5 and group 6 are experimental hornblende peridotite 32 and garnet pyroxenite 
melts 36,37 , respectively. Olivine fractionation and accumulation, 1 wt% step. Fractional (FC) and equilibrium (EC) clinopyroxene crystallization, 10 wt% 
step. The FC3MS values of the melts of the three source rocks are estimated by the relationships between FC3MS value and melting degree shown in 
Fig. 2c. The La/Yb ratios are caculated by batch melting model, modeling parameters see Supplementary Table 6. The correlation coefficients (R 2 ) 
calculated using logarithmic line. The assumed error bars of the natural basalts represent ± 10% variation of the FC3MS values and La/Yb ratios, which 
probably result from many secondary factors including partial melting conditions, magma evolution, analytical error and source compositional diversity 
etc. Other explanations same as in Fig. 1. 



similar as fractionation trends of olivine and clinopyroxene and 
plagioclase, however, which can also be attributed to the contribu- 
tions from the pyroxenite melts (Fig. 3a), but currently it is difficult to 
verify and estimate the proportion of the pyroxenite melts because of 
the wide variations in the melts compositions for the heterogeneous 
peridotite and pyroxenite sources. 

In Fig. 3b, for peridotite model, the negative correlations between 
the FC3MS value and MgO content of the selected four groups C- 
OIB 273031 can be readily explained as the result of different degrees 
fractionation of olivine and clinopyroxene, and the group 3 and 



group 4 data seem result from clinopyroxene-dominated fractiona- 
tion. However, pyroxenite melts (group 6) also can produce these 
trends. And, it is difficult for such high MgO melts (for example, 
group 3, MgO = 8.6-10.8 wt%; group 4, MgO = 11.3-12.8 wt%) to 
crystallize clinopyroxene. The most important evidence to support 
pyroxenite source is that the significant positive correlations between 
the FC3MS value and La/Yb ratio for each group data are in good 
agreement with the melting trend of model average pyroxenite 
(Fig. 3c). The pyroxenite melting trends can only be explained 
by crystallization process of peridotite melts when assuming that 
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clinopyrxene or other fractionating minerals (for example, plagio- 
clase) have extremely high mineral/melt partition coefficient of Yb, 
just as garnet, but this is rather unrealistic. In addition, the horn- 
blende peridotite, may represent metasomatized peridotite, cannot 
produce high FC3MS value melts, although it can produce high La/ 
Yb ratio melts (group 5) 32 . The simultaneous consideration of 
FC3MS versus MgO and FC3MS versus La/Yb in Fig. 3b, c suggests 
that many reputed evolved low magnesian (MgO, —4-10 wt%, Mg # , 
~50-60) C-OIBs in fact represent primary pyroxenite melts. These 
low magnesian primary magmas might also be present in various 
types of basalts such as typical ocean island basalt and continental 
flood basalt, although they are commonly deemed as evolved mag- 
mas. Thus it should be cautious about using olivine-correction 
method to gain high magnesian primary magmas that equilibrium 
with typical mantle peridotite. One important consequence is that 
corrected high magnesian magma would overestimate mantle poten- 
tial temperature, a similar conclusion had been reached by Presnall 
and Gudfinnsson 17 . 

Phase relation constraints on the C-OIB. Trace elements suggest 
that the C-OIB have the typical garnet signature, which commonly 
are attributed to various degrees mixing between garnet Iherzolite 
melts and spinel Iherzolite melts 25,29 . This explanation seems to be 
supported by projections in phase diagrams where the C-OIB, except 
for some highly alkaline basalts, overlap with volatile-free peridotite 
melts (Fig. 4). Note that olivine fractionation in Fig. 4a will increase 
the overlap between the C-OIB and volatile-free peridotite melts. 
Nevertheless, all the peridotite melts cannot explain the high 
FC3MS values and the relationship between the FC3MS values 
and La/Yb ratios of the C-OIB (Fig. 3). Therefore, the projections 
of the C-OIB in Fig. 4 indicate that they are mostly partial melts 
of garnet pyroxenite formed on the cotectic [L + Cpx + Gt]. Note 
that the cotectic [L + Cpx + Gt] may also contain small amount 
of spinel, olivine and orthopyroxene, because of variations in 
source bulk composition, melting degree and pressure. Note also 
that the 1.6-2 Gpa cotectics [L + Ol + Cpx + Opx + Gt] of 
hydrous peridotite are similar as the alkaline C-OIB projections, 
whereas the 3 Gpa cotectic [L + 01 + Cpx + Gt] of hydrous 
peridotite richer in diopside component than that of volatile-free 
peridotite (Supplementary Fig. 26b). Although hydrous peridotite 
at 1.6—2 Gpa can make garnet stable and generate melts that 
overlap with the C-OIB projections in the phase diagrams (Fig. 4), 
the FC3MS values of these hydrous melts are too low (- — 0.12- 
0.01) 33 to explain the high FC3MS values of the C-OIB. Therefore, 
hydrous peridotite should not be a good candidate for the C-OIB 
source. 

For the garnet pyroxenite melts formed on the cotectic [L + Cpx 
+ Gt], one important characteristic is that the melts projections no 
longer record their detailed generation pressures (Supplementary 
Fig. 27). The tholeiitic C-OIB might be generated at similar or even 
higher pressures than the alkaline C-OIB. These pressure informa- 
tion acquired from the pyroxenite model are significantly different 
from those based on peridotite model, where alkaline basalts are 
typically assumed to be generated at higher pressures and lower 
melting degrees than tholeiitic basalts 24,25,2 *. In addition, even for 
the cotectic [L + Ol + Cpx + Gt ± Opx], the pressure information 
still can be significantly affected by addition of C0 2 , but these melts 
have too much diopside component than the C-OIB (Supplementary 
Fig. 28). The highly silica undersaturated C-OIB are also very similar 
as carbonated garnet pyroxenite- and hornblendite-derived melts 
(Supplementary Figs. 27 and 28), indicating that their source might 
contain C0 2 , hornblende and phlogopite etc. 

Discussion 

The best and our preferred explanation is that the C-OIB are mainly, 
if not all, pyroxenite melts. The pyroxenite melts with low Mg* and 
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Figure 4 | Projections (mole%) of the C-OIB compared with 
experimental melts that equilibrium with garnet pyroxenite and 
peridotite. (a), A projection from or towards Diopside into the plane 
Olivine-Quartz-Calcium Tschermak's. (b), A projection of the same data 
as in (a) from or towards Olivine into a part of the garnet-pyroxene plane. 
This pyroxene-garnet plane is represented as the Opx-CATS line in (a). 
Projection code is from ref. 2 (see Supplementary Fig. 4). The regions 
depicted by dashed line are the potential basaltic melts that equilibrium 
with garnet pyroxenite, and the red arrows indicates the direction of 
decreasing temperature of melts that equilibrium with major residual 
phases, garnet and clinopyroxene. The peridotites shown in (a) and (b) are 
the source rocks of the volatile-free peridotite melts (>2 Gpa). The 
detailed projections presented in Supplementary Figs. 4, 26 and 27. CATS- 
CaAl 2 Si0 6 , Diopside- Ca(Mg, Fe)Si 2 0 6 , Olivine- (Mg, Fe) 2 Si0 4 , Opx 
(Enstatite)- (Mg, Fe) 2 Si 2 0 6 , Quartz- Si0 2 , Grossular- Ca 3 Al 2 Si 3 0 12 , 
Pyrope- (Mg,Fe) 3 Al 2 Si 3 0 12 , Wo- CaSi0 3 , Flp- (Na, Ca)(Si, A1) 4 0 8 or (Na, 
K)AlSi 3 0 8 , Rut-Ti0 2 , Sp-MgAl 2 0 4 . 

high alkali contents are experimentally characterized by very low 
temperatures, even lower than the normal mantle peridotite solidus 6 . 
When the pyroxenite partially melt to produce basaltic magmas that 
have similar compositions as the C-OIB, the mantle peridotite may 
still subsolidus, the pyroxenite melts, however, can have higher melt- 
ing degree than commonly assumed peridotite melts (Supplemen- 
tary Information), and thus might have much stronger migration 
ability. This probably is one of the most important reasons that the 
C-OIB, without equilibrium with mantle peridotite, can record the 
characteristics of pyroxenite melts although carrying a variety of 
crustal and mantle xenoliths. Therefore the garnet pyroxenite would 
partially melt when mantle decompression or asthenosphere upwel- 
ling occurs at normal mantle potential temperature. This process 
does not need the high temperature mantle plume model 34 (Supple- 
mentary Information). 
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The C-OIB have broadly chemical similarities and small in volume 
although temporally and spatially widely distributed 24,25,29 . This is 
likely to be the characteristic features of continental pyroxenite- 
derived melts as a result of mantle lithological heterogeneity after 
long periods recycling of continental crust 8 . The two-component 
marble-cake mantle model 4 or the statistical upper mantle assem- 
blage model 7 might better characterize the continental upper mantle 
and seems to be supported by the pyroxenite-derived C-OIB. 
Importantly, because there are significant differences with regard 
to the petrogenesis of peridotite- and pyroxenite-derived basalts, 
many previous geological and petrological interprations of the C- 
OIB based on the peridotite model 24,25,29 need to be reconsidered. 

In conclusion, we propose that the FC3MS value is a useful cri- 
terion for identification of source lithology of basalts even if taking 
into account source compositional diversity and variation of melting 
and crystallization conditions of basaltic magmas. This discriminant 
parameter provides a new perspective to understand the composi- 
tional characteristics of basalts and also provides insight into the 
origin of basalts. 

Methods 

In trying to find the taxonomical criterion of source lithology of basalts, we have 
collected nearly all the experimental melts on a variety of peridotite and pyroxenite in 
literature before the year 2013. Note that our selected peridotite melts are not 
exclusively on volatile-free or normal mantle peridotite. All the experimental melts 
and C-OIB are carefully examined using simple plots, projections and FC3MS-based 
diagrams. The parameterization results of several representative perdotite and pyr- 
oxenite melting trends are presented and explained in Supplementary Figures 8-11. 
Other details about the data selection and sources of experiments and the C-OIB and 
the parameterization are presented in Supplementary Information and Figures 1-3. 

We do not assume either peridotite or pyroxenite end-member source composition 
for mantle- derived magmas, but they can be roughly divided into two lithological 
units using their whole-rock MgO contents. Although the majority of the experi- 
ments are batch melting results, our parameterized FC3MS value would not change 
significantly for accumulated fractional melting (Supplementary Fig. 22). 
Importantly, our main conclusion that garnet pyroxenite is the dominant source 
lithology for the C-OIB is based on a safely upper limit of FC3MS value —0.65, which 
approximate to near-solidus melts of iron-rich peridotite. If peridotite melts are 
mixed products from a variety of peridotite under various conditions (e.g. volatile- 
free and volatile present), the resultant melts still have FC3MS values less than 0.65 
(Supplementary Information). Nevertheless, the magma mixing mechanism would 
easily change the major element content, which might be the reason that peridotite 
and pyroxenite melts cannot be distinguished in the commonly used Harker diagrams 
and phase diagrams. Although olivine chemistry method 1019 has been widely used to 
identify pyroxenite melts, it might introduce many uncertainties due to source 
compositional diversity and melting and crystallizing conditions before enough high 
precision experiments achieved (Supplementary Information). For some low-mag- 
nesium C-OIB, high FC3MS, high Fe/Mn ratios and low CaO contents might result 
from fractionation of clinopyroxene and olivine, but pyroxenite melting trends can be 
indicated by FC3MS and La/Yb relations (Fig. 3). 
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